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Abstract. The topological structure of the QCD vacuum can be probed by monitoring the spatial 
localization of the low-lying Dirac eigenmodes. This approach can be pursued on the lattice, and 
unlike the traditional one requires no smoothing of the gauge field. I review recent lattice studies, 
attempting to extract a consistent description. What emerges is a picture of the vacuum as a 
"topological sandwich" of alternating, infinitely thin 3d layers of opposite topological charge, as 
originally seen in direct measurements of the topological charge density. 
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"Understanding" confinement, by identifying the relevant infrared degrees of freedom 
of the gauge field, has been a long-standing theoretical goal. It is natural to associate 
this non-perturbative phenomenon with non-perturbative, topological excitations. The 
standard list of potentially relevant excitations consists of instantons, Abelian monopoles 
and center vortices, with co-dimension 4, 3, and 2 respectively. They each have received 
a fluctuating degree of attention over the years. One may hope that a proper lattice 
study may unambiguously identify the right excitation. However, in the past it has been 
necessary to filter out UV fluctuations of the gauge field in order to reveal the large- 
scale structure. This is accomplished by a smoothing/cooling/smearing procedure which 
reduces the action, and inevitably drives the gauge field towards an action minimum, 
i.e. a classical instanton solution. Recently, as outlined below, a different strategy has 
been followed, which avoids such bias. The localization properties of low-lying Dirac 
eigenmodes presumably tell us about the underlying gauge field excitations, responsible 
for chiral symmetry breaking and confinement. 

LOCALIZATION: ANDERSON AND DIAKONOV-PETROV 

Anderson [Q]] considered the Hamiltonian H = A + V, where A is a nearest-neighbor 
hopping operator (a discretized Laplacian) and V a random potential. This Hamiltonian 
mimics that of a crystal doped with random impurities. An eigenmode \jf(f) can be 
localized, meaning that | yr(r) | decays exponentially for large |r| : the electron cannot hop 
to infinity, and this mode does not contribute to the electric conductivity of the material. 
Otherwise, the eigenmode is extended. Anderson showed that eigenmodes were always 
localized, if the disorder in V was large enough, or the energy low enough. This is 
intuitively clear: at very high energy, the random potential plays no role and eigenmodes 
are extended, plane-wave-like; but at low energy, sufficient randomness in the potential 
may forbid hopping to any of the neighboring sites. Thus, the spectrum looks genetically 
as in Fig. 1, with a mobility edge X c separating the localized from the extended regime. 
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FIGURE 1. Localization properties of boson eigenstates in a random medium. 



X c may lie above or below the ground-state energy depending on the disorder, allowing 
for a so-called quantum transition at zero temperature. 

Diakonov and Petrov [2] proposed an analogous explanation for the QCD chiral tran- 
sition at temperature T c . Since (i/ai/a) = — % p(0) 030, the spectral properties of the Dirac 
operator, described by the density p(A), must change with temperature. Recalling that 
an instanton supports a chiral zero mode i/a 7 (for each fermion flavor), one considers the 
Dirac spectrum for a linear superposition of Instantons and Antiinstantons. The exact 
zero modes are now displaced. The infrared spectrum results from diagonalizing the ef- 



ual Instantons and Antiinstantons zeromodes. At T = 0, this overlap decreases as 



and gives rise to extended modes having support on all instantons and antiinstantons, and 
an essentially uniform spectrum near zero. Such extended modes can be observed on the 
lattice (see Fig. 2). As the temperature increases, the density of instantons decreases 
(their action 87T 2 /g 2 (T) increases) and now decays exponentially ~ exp(— Ttr^T) in 
spatial directions. Both factors may trigger a transition to localization, which suppresses 
near-zero eigenvalues and restores chiral symmetry. 

As in Anderson's bosonic case, this is a disorder-driven transition. But Dirac eigenval- 
ues come in pure imaginary pairs =b'A, and the focus here is on the center of the spectrum 
A = 0. Correspondingly, the transition can be modeled by chiral random matrices [50. 

It is important to note that the Diakonov-Petrov scenario does not require instantons, 
but only chiral zero modes. It turns out that other topological defects - domain-walls, 
monopoles, vortices - also support chiral zero modes. Under the working assumption 
that extended modes have support on the union of topological defects, this opens the 
possibility to determine the topological vacuum structure from the spatial distribution 
of low-lying eigenmodes. This approach is gauge-invariant and requires no smooth- 
ing/cooling of the gauge field. 



To characterize how localized or extended a mode \\f{x) is, one uses the inverse partici- 
pation ratio I PR = V (^i^^w - This ratio of moments is equal to 1 if yr(x) = 8(xq) is 

completely localized, and to V if i/a(x) = const. Wx is completely delocalized. If i/a(jc) = 1 
on a fraction / of the sites, elsewhere, then I PR = 1/ /, which justifies its name. 
Of course one should consider the continuum limit a — > of the lattice study. Note 




LATTICE STUDIES 




FIGURE 2. Surface (yellow) of equal magnitude of lowest-lying Dirac eigenmode in the presence of 
an Instanton-Antiinstanton pair (red and blue spheroids), from 0]. The eigenmode is extended over both 
objects. The gauge field was cooled to identify the I — A pair, but not to obtain the Dirac eigenmode. 




FIGURE 3. The fraction / of occupied lattice sites, hence the IPR, remains unchanged as a — ► in the 
presence of any thick object (left). It goes to zero as a d ~ 4 for a thin object of dimension d (right). 



that the fraction / of occupied sites, hence the IPR, remains constant as a — ► for 
any type of macroscopic object. Conversely, if the object is "thin" and lives on a 



submanifold of dimension d and volume men / = yy^t ~ 07 "> ana Lt " K ~ a " 
which diverges as a — > (see Fig. 3). d = 0, 1,2 then characterizes "thin" instantons 
(point-like), monopoles (line-like) and vortices (surface-like). 

• The first study, by the MILC collaboration [jg] for SU(3), indicated d ~ 3: the 
vacuum, it seems, is made of infinitely thin domain- walls! (see Fig. 4, left). 
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FIGURE 4. (left) Most recent IPR data from MILC Collaboration [6] versus lattice spacing, for SI/ (3) 
with Asqtad Dirac operator and Symanzik gauge action. The dimension of the supporting manifold is 
about 3. (right) Correlator (| l/f(0)||v(x)|) of magnitude of low-lying Dirac eigenmode, for different lattice 
spacings. The straight line shows what would happen if = 1 on a 3d fractal, elsewhere. 




• This surprising result was quickly checked in [7]. There, the IPR diverges even 
faster, favoring point-like instantons (see Fig. 5). That study used the overlap discretiza- 
tion of the Dirac operator, with exact chiral properties, but used the Wilson action for the 
gauge group SU (2) . Unfortunately, it is known since Q8D that with this choice of gauge 
action, the action of size a defects ("dislocations") is insufficient to compensate for their 
entropy ~ loga~ 4 , and they become dense as a — > 0. The observations of 10] may then 
be entirely caused by lattice artifacts. 

• A third study, also with the overlap Dirac operator, but for SU (3) and Liischer- 
Weisz gauge action, shows results (Fig. 6, left) roughly consistent with [6]. The IPR for 
the low-lying modes grows more or less as I /a, consistent again with thin domain- wall 
structures (Fig. 6, right). 

It should not have come as a surprise that the IPR diverges. Even if classical lumps of 
size I/Aqcd are relevant to the QCD vacuum structure, these lumps do not look smooth 
at short distance due to quantum fluctuations (see Fig. 7). Besides, this result can be 
related to a curious property of the topological charge density operator (q(0)q(x)) fl 1 ON - 
On one hand, this correlator is negative for any x ^ 0. This can be seen from reflection- 
positivity, or simply by realizing that q(x) ~ E ■ B acquires an extra factor i in Euclidean 
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(left) IPR versus eigenvalue and lattice spacing, for SU(3) with overlap Dirac operator 
Weisz gauge action [9]. (right) For the lowest non-zero modes, the IPR grows oc \ja as for 




FIGURE 7. A classical kink is smooth. The actual kink looks like a random walk at short distances. 



space. On the other hand, (/ d 4 x q(0)q(x)} = Xtop is finite and about (190MeV) 4 . 
Hence, a contact term 5 (3c) must compensate the negative x ^ integral. Moreover, 
the canonical dimension of q(x) is 4, so that one expects (J^ d 4 x q(0)q(x)) to diverge. 
This divergence is exactly compensated by the contact term to leave the finite piece Xtop- 
This behaviour has recently been observed on the lattice: compare Fig. 8 left and right. 



Topological charge density correlator 
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FIGURE 8. Correlator of the topological charge density: (left) as expected in the continuum and (right) 
as measured on the lattice [11]. The negative part and divergent contact term both become visible as a —* 0. 



FIGURE 9. Topological charge density: (left) in QCD H and (right) in the (1 + \)d CP 3 model fill . 



One has then to understand how the vacuum becomes "topologically antiferromagnetic", 
with (q(0)q(x)) going to — °° as |x| — > 0. Before the studies of Dirac eigenmodes 
reviewed here, the Kentucky group already claimed to see space-filling 3d structures 
of transverse size ff[a) and of opposite topological charge s 111 ill . Analogous structures 
have now been seen in Ref. D90, and also in the CP 3 model I112n (see Fig. 9). 

As a final check, the MILC collaboration has measured the correlator ( | y (0) 1 1 y(x) | ) , 
where y is a low-lying Dirac eigenmode. If | ty(x) | = 1 on a 3 d fractal and elsewhere, 
this correlator would decay as l/\x\. Measurements are not too different (Fig. 4, 
right), until |3c| reaches distances Kqqd). 

• Most recently, this kind of lattice study has been extended to bosonic fields and 
eigenmodes of the co variant Laplacian in various representations 1130, with yet more 
surprises. In the adjoint representation, the mobility edge is found to rise unexpectedly 
as 1/a (see Fig. 10). If true, all finite-energy modes in the continuum would be localized, 
contradicting perturbation theory. One wonders if this strange result also can be blamed 
on the dense dislocations caused by the choice of Wilson action for gauge group SU (2) . 
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FIGURE 10. The mobility edge (arrow) for an adjoint Higgs appears to go to °° as a — > 11311 . for 
SI/ (2) with Wilson gauge action. 




FIGURE 11. The "topological sandwich" structure observed at short distance (blue) is not inconsistent 
with another structure (red) at scale 1 /Aqcd. 



CONCLUSION 



The picture of the QCD vacuum which emerges (Figs. 4, 6, 8, 9) is that of a "topological 
sandwich", with alternating, infinitely thin 3d layers of opposite topological charge 
density. While bizarre, this picture is not forbidden, and is supported by some theoretical 
arguments [10]. But it is far from the usual picture of a dilute gas of classical excitations. 
One should stress that these two descriptions are not mutually exclusive, as sketched 
Fig. 11: one applies at UV scales, the other may apply at scales &{\/ Aqcd)- The 
structure of the QCD vacuum depends on the scale considered. 
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